Virial expansion for a strongly correlated Fermi gas with imbalanced spin populations 
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Quantum virial expansion provides an ideal tool to investigate the high-temperature properties of 
a strongly correlated Fermi gas. Here, we construct the virial expansion in the presence of spin pop- 
ulation imbalance. Up to the third order, we calculate the high-temperature free energy of a unitary 
Fermi gas as a function of spin imbalance, with infinitely large, attractive or repulsive interactions. 
In the latter repulsive case, we show that there is no itinerant ferromagnetism when quantum virial 
expansion is applicable. We therefore estimate an upper bound for the ferromagnetic transition 
temperature T c . For a harmonically trapped Fermi gas at unitarity, we find that (T c ) uppper < Tf, 
where Tf is the Fermi temperature at the center of the trap. Our result for the high-temperature 
equations of state may confront future high-precision thermodynamic measurements. 



I. INTRODUCTION 

Strongly interacting Fermi gases of ultracold atoms 
near Feshbach resonances have offered unique opportu- 
nity to study some long-sought problems in condensed 
matter physics, astrophysics and particle physics [HE!- 
An particular example that attracts intense interests is 
the ground state of a two-component, spin-population 
imbalanced Fermi gas of 6 Li or 40 K atoms [!, H|. In 
the strongly interacting regime, observation of new ex- 
otic superfiuid states in the imbalanced systems may be 
anticipated, such as the spatially inhomogeneous Fulde- 
Ferrell-Larkin-Ovchinnikov state 0-[T2|- However, be- 
cause of strong interactions, theoretical investigation of 
imbalanced Fermi systems turns out to be difficult and 
restrictive. 

In this work, we theoretically investigate the high- 
temperature equation of state of a strongly interacting, 
imbalanced Fermi gas, by generalizing a quantum virial 
expansion method to imbalanced systems (l3l - [l6l . lill [l9| . 
Virial expansion for balanced Fermi gases with equal spin 
populations has already shown to be very useful, in the 
studies of both static p~3MT7l | and dynamic properties 
fl8L f]~^ | . The third virial coefficient of a unitary Fermi gas 
with infinitely large scattering length has been calculated 
[l5| and consequently confirmed experimentally (2^|. The 
third virial coefficient of strongly interacting Bose gases 
has been determined [21]. The virial expansion of single- 
particle spectral function [l8[ or dynamic structure fac- 
tor [l9[ has also been developed to understand, respec- 
tively, the relevant measurements of momentum-resolved 
rf spectroscopy (22l. [23| or Bragg spectroscopy (24J. 

Our generalization of quantum virial expansion to the 
imbalanced systems gives an approximate means to study 
interesting phenomena associated with the spin degree 
of freedom. For instance, for a strongly repulsively in- 
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teracting Fermi gas, itinerant ferromagnetic transition is 
expected at low temperatures (25l - [30l |. However, the crit- 
ical transition temperature is yet to be determined. The 
use of virial expansion may allow us to estimate an up- 
per bound for the critical temperature. By expanding 
to the third order, we find that there is no itinerant fer- 
romagnetism as far as the quantum virial expansion is 
applicable. Therefore, we show that for a strongly repul- 
sive Fermi gas in harmonic traps, the critical temperature 
should be smaller than the Fermi temperature at the cen- 
ter of the trap. We note that, virial expansion method 
has been used by Horowitz and Schwenk to study the low 
density, spin-imbalanced neutron matter [3l[ l32y. Their 
expansion was taken up to the second order, where the 
second virial coefficient was calculated from the nucleon- 
nucleon scattering phase shifts. 

Our paper is organized as follows. In the next section, 
we show how to expand the thermodynamic potential 
(or pressure in case of a uniform system) of an imbal- 
anced Fermi system in terms of virial coefficients. The 
n-th virial coefficients can be calculated from the parti- 
tion functions of a cluster that contains up to n-particles. 
Different spin configurations in a n-particle cluster give 
rise to different virial coefficients. In Sec. Ill, we focus on 
the expansion to the third order and calculate the finite 
temperature free energy of an imbalanced Fermi gas as 
a function of spin imbalance. The possibility of itinerant 
ferromagnetism in a strongly repulsively interacting sys- 
tem is discussed. Sec. V is devoted to conclusions and 
some final remarks. 



II. QUANTUM VIRIAL EXPANSION FOR AN 
IMBALANCED FERMI GAS 

The starting point of quantum virial expansion is that 
at large temperatures the chemical potential /i diverges 
to —oo. Therefore, the fugacity z = exp(fi/ksT) = 
exp(/3//) <C 1 is a well-defined small expansion parameter. 
We may expand the high-temperature thermodynamic 
potential Q of a quantum system in powers of fugacity, 
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no matter how strong the interaction strength would be. 

For a two-component Fermi ga s with balanced spin 
populations, it was shown that (la. Ilijj]. 

= -k B TQ 1 [z + b 2 z 2 + ■■■+ b n z n + ■■■], (1) 

where Q n — Tr„[exp(— "H/fcsT)] is the partition function 
of a cluster containing n particles and b n is the n-th virial 
expansion coefficient. The trace Tr n takes into account 
all the spin-configurations of n-particles and traces over 
all the states of a proper symmetry. The virial coefficient 
b n was found to take the form (la . EH , 



Qa/Qi 



l/2, 



X2 



and 



x 3 = z^Q 3fi + z*ZiQ 3 ,i + z t zfQ 3 . 2 + zfQ3,3- (5) 



By introducing a symmetric cluster partition function 
Q n = Q n<0 = Q n}n and using the properties of Q„ t k, 
it is easy to show that Q\ — Qi/2, Q 2 ,i — Q2 — 2Q 2 , an d 
<?3,i = = Q3/2 - Q|. Using ln(l + x\ + x 2 + X3) ^ 
(xi +X2+X3) — (x\ + 2xix 2 ) / 2 + x\ /3, after some algebra 
we obtain b n ,k {k < n/2), 



and 



61,0 = 1/2, 

^2,0 = Q2I Qi ~ Qti 



(G) 
(7) 



&3 = Q 3 /Qi-Q 2 +Qi/3, etc. 



(2) 



In the presence of spin imbalance, it is necessary 
to introduce two fugacities Zf = exp^-f/fc^T) and 
zj. = exp(/i^/fcsT'), and to distinguish different spin- 
configurations. Quite generally, we may write the ther- 
modynamic potential as, 



n=l k=Q 



~ k 7 k h 
Z 1 O r , 



(3) 



where 6 n ,fe is the n-th (imbalanced) virial coefficient con- 
tributed by the configuration with n — k spin- up fermions 
and k spin-down fermions. It is easy to see that the 
imbalanced virial coefficients satisfy the relation b n ^ — 

and YJl=Q b ™,k = 
The calculation of b n .k is straightforward, following 
the standard definition of thermodynamic potential. We 
rewrite the grand partition function Z = TYexp[— (Ji — 
\x^N^ - h^M^/UbT] in the form, 



2 = E£*T**ftr 



(4) 



n=0 fe=0 



where is the partition function of a cluster that con- 
tains n — k spin-up fermions and k spin-down fermions. 
It is apparent that due to the symmetry in spin con- 
figurations we have Q n .k — Q n ,n-k- The imbalanced 
cluster partition functions satisfy as well a sum rule 
Sfe=o Qn.k — Qn- By expanding the thermodynamic po- 
tential £1 = — ksTlnZ into powers of the two fugacities, 
the imbalanced virial coefficients can then be expressed 
in terms of the cluster partition function Q n ^,- 



A. Virial expansion up to the third order 

To be concrete, let us consider the imbalanced virial 
expansion up to the third order. To this order, we may 
write the grand partition function as Z = l+xi+X2 + x 3l 
where 



xi = ZfQi t o + z±Qi,i, 



62,1 = Q2/Q1 ~ 2Q S 2 /Qi - Qi/4, 



Q3M ~ Q2/2 + Qi/24, 



(8) 
(9) 



and 



J3,i 



Q 3 /(2Qi 



I/Q1-Q2/2 + QI/2 + QI/8. (10) 



The virial coefficients with k > n/2 can be obtained di- 
rectly since b n . k = b n . n - k . 

In practice, it is convenient to consider the interac- 
tion effect on the virial coefficients or the differences 
such as AQ n = Q n - Ab n = b n - b£ , and 

Ab n ,k — b nt k — b^\. Here, the superscript "1" denotes 
an ideal, non-interacting system with the same fugacities 
and the operator "A" removes the non- interacting con- 
tribution. It is clear that the symmetric cluster partition 
function Q s n is not affected by interactions since the in- 
teratomic interaction occurs only between fermions with 
unlike spins. Thus, we have Afe2,o = A^o — 0, A&2,i = 
A(6 2 -26 2 ,o) = A6 2 and A6 3 ,i = A(6 3 /2-6 2 ,o) = A6 3 /2. 
Accordingly, we may rewrite the thermodynamic poten- 
tial into the form (up to the third order), 



n = - ksTQi 



z- t z i Ab 2 H Ab 3 



, (11) 



where f^ 1 ' is the thermodynamic potential of a non- 
interacting Fermi gas. 



B. Thermodynamic potential Si' 1 ' 

The thermodynamic potential of an ideal, imbalanced 
Fermi gas is simply the sum of thermodynamic poten- 
tial of each component, flW — fiW^if) + Q^(u^). In 
homogeneous space, where Qi : h = 2V/X 3 , the single- 
component thermodynamic potential 0^'(/x ff ) takes the 
form (a =t,l), 



n« = -v kBT 2 



A 3 



i 1/2 In (l + z^e"*) di, (12) 
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where V is the volume and A = [2nh 2 / (mfcflT)] 1 ^ 2 is the 
thermal wavelength. We have used the subscript U H" to 
denote the result in the homogeneous space. 

In the presence of a harmonic trap Vr(r) = mwV/2, 
we consider the thermodynamic limit with a large num- 
ber of fermions. The non-interacting thermodynamic po- 
tential f2^P(/v) is then given semiclassically by (/3 = 
l/k B T), 



1 f drdk 



In 



l + e 



PJ (2tt)' 

^) I / /-In (1.4 ;. T( ')'//• 
(M 2 



(14) 



Accordingly, we have Qi,t = 2 (ksT) 3 / (hui) 3 . Here, the 
subscript "T" denotes the result in a harmonic trap. 



C. Virial coefficients A62 and Ab3 

The calculation of virial coefficients requires the knowl- 
edge of full energy spectrum of a few-body system. By 
utilizing the exact two- and three-fermion solutions, we 
have calculated the second and third virial coefficients for 
a strongly attractively or repulsively interacting Fermi 
gas [16j . Here, the terminology of "a strongly attrac- 
tively interacting Fermi gas" refers to a resonant Fermi 
gas near a Feshbach resonance with the inclusion of the 
molecular states, while the terminology of "a strongly re- 
pulsively inteacting Fermi gas" means a resonant Fermi 
gas with the molecular branch excluded. In other words, 
all fermions in the strongly repulsively interacting Fermi 
gas populate only on the upper, atomic branch of the 
energy levels near a Feshbach resonance. At unitarity, 
three resonantly interacting fermions in an isotropic har- 
monic trap can be solved exactly using hyperspherical 
coordinates [33|, |34j, with the relative energy given by 
E re i = (2q + si t n + l)fku, where q and si_ n (I the rela- 
tive angular momentum and n the number of nodes) are 
respectively the good quantum numbers in the hyperra- 
dius and hyperangle equations. One can show that the 
molecular branch corresponds to the states with s^ n= Q 
[lfjj . These states should be excluded for the "repulsive" 
case. We refer to the Sec. IIIB of Ref. [l6[ for a detailed 
discussion about this. 

In the following, we shall focus on the unitarity limit 
with infinitely large scattering length, where the virial 
coefficients are universal and temperature independent. 
The coefficients for a homogeneous gas are given by [l6[ , 



Auatt _ 
^°2,H — 



+1/V2, 
-1/V2, 



(15) 
(16) 



Here, the superscript "att" (or "rep") means the coeffi- 
cient of an attractively (or repulsively) interacting Fermi 
gas. The third virial coefficient of an attractive unitary 
Fermi gas, Ab 3 u H ~ —0.355, was recently confirmed ex- 
perimentally at Ecole Normale Superieure (ENS), Paris 

For a trapped Fermi gas, we instead have, 



A6g£ = 



4/4, 
-1/4, 



and 



AblJ 



-0.06833960, 
0.34976. 



(19) 
(20) 

(21) 
(22) 



We note that, in the unitarity limit the homogeneous 
and trapped virial coefficients are connected universally 
by the relation (TH llfjj. 



.H 



,3/2- 



(23) 



The factor of n~ 3 / 2 reduction in harmonic traps is sim- 
ply due to the higher density of states in traps. To see 
this, let us consider the thermodynamic potential of a 
harmonic trapped Fermi gas in the local density approx- 
imation f2 = J drfi(r), where fl(r) is the local thermo- 
dynamic potential 



fi(r) oc z (r) + b 2M z 2 (r) + b 3M z 3 (r) + 



(24) 



and 



Here, the local fugacity z (r) = zexp[—V(r)/kBT] is de- 
termined by the local chemical potential /i(r) = fi — V(r). 
The spatial integration immediately leads to Eq. (f2"3"| . It 
is obivious that the suppressed virial coefficients in har- 
monic traps imply a weaker interaction effect and conse- 
quently a better convergence for virial expansion. 



D. Virial expansion of high-T free energy 

We are now ready to study the high temperature 
thermodynamics of an imbalanced, strongly interacting 
Fermi gas. With the virial expansion of thermodynamic 
potential Eq. (jlll) , we solve the standard thermodynamic 
relations iVf = —dfl/dfpf and iVj, = —dVL/dfjL^ for the 
two fugacities z-j- and z±, at a given reduced temperature 
t = T/Tp and a given spin imbalance P = (Nf — Ni)/N. 
Here, Tp is the Fermi temperature. It is given by 
Tp = ft 2 (37r 2 n) 2 / 3 /(2m)/kB in the homogeneous space 
and by Tp = {ZNf^hu/kB in a harmonic trap. 

In the homogeneous space, it is easy to show that we 
can define a dimensionless number density h — n\ 3 = 
8/(3v^Ft 3 / 2 ), hf = (1 + p)n/2, and n x = (1 - p)n/2. 
We then rewrite the number equations into dimensionless 
forms, 



Abl% 



-0.35501298, 
-1.8174. 



(17) 
(18) 



n t = n (z t ) + z t z i 2Ab 2 + (2z 2 z i + z^zf) Ab 3 , (25) 
H = n (zi) + z tZl 2Ab 2 + (z^z i + 2z t zl) Ab 3 , (26) 
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where n(z) = (2/ y /Hr)f™rf[ze- t /(l + ze- t )]dt. We 
can obtain the two fugacities by solving the coupled 
number equations. Then, we calculate the free energy 
F = ft + [4fNf + HiN^. The free energy per particle in 
units of Fermi energy Ep is given by, 



F 



NEt, 



l+p 



In 



1-p 



In zi 



2 1 2 
+/ {h) + ^z i 2Ab 2 + (z^zi + z t zf) Ab$,7) 

where / (z) = (2/y / 7r) J °° y/thi (1 + ze~ l ) dt and the con- 
stant A = 3V7tt 5 / 2 /8- 

In a harmonic trap, we have very similar equations 
with exactly the same structure. However, the dimen- 
sionless number density is now given by, n — l/(3r 3 ). 
The function n(z) and f{z) shall take the form, 



n (z) 



-t'dt, 



(1 + ze~*) 
1 C°° 

f(z) = - J \n(l + ze- t )t 2 dt, 



(28) 
(29) 



respectively. Finally, the constant A in free energy should 
be read as A = 3r 4 . 



III. RESULTS AND DISCUSSIONS 

We discuss separately the high-temperature free energy 
of a unitary Fermi gas in homogeneous space and in a 
harmonic trap. 



A. Free energy of a homogeneous unitary Fermi gas 

Fig. 1 presents the virial expansion prediction for the 
free energy as a function of temperature, for a balanced 
unitary Fermi gas or for a spin imbalanced system. For 
comparison, we show also the free energy of an ideal, 
non-interacting Fermi gas, .To illustrate clearly the 
interaction effect, we show in Figs, lb and lc the residue 
free energy or the interaction energy AF = F - fW, 
by subtracting the non-interacting background F^K As 
shown in Fig. la and in the lower part of Fig. lb for a 
balanced Fermi gas, in the presence of attractive interac- 
tions the second- and third-order virial expansion predict 
quantitatively the same free energy down to the Fermi 
degenerate temperature Tp. This suggests a very broad 
temperature window for applying the quantum virial ex- 
pansion method. In case of repulsive interactions, how- 
ever, the good agreement from the expansions at differ- 
ent order is found only at T > 5 If, indicating that the 
virial expansion for a repulsively interacting Fermi gas is 
severely restricted. The similar applicability of virial ex- 
pansion is observed as well for an imbalanced Fermi gas 
with a general spin imbalance P. On the other hand, by 
comparing the interaction energy at P = (Fig. lb) and 
P = 0.5 (Fig. lc), a finite spin imbalance reduces the 
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Figure 1: (color online) Temperature dependence of the free 
energy of a homogeneous unitary Fermi gas, calculated up 
to the second-order (dashed lines) or third-order (solid lines) 
using the quantum virial expansion theory. For comparison, 
we show the ideal, non-interacting free energy by a thin dot- 
dashed line. As anticipated, the free energy of an attractive 
or a repulsive Fermi gas is always smaller or larger than that 
of an ideal Fermi gas, respectively. 



interaction energy as expected, although the reduction is 
not so significant. 

Fig. 2 reports the free energy as a function of the spin 
imbalance at a given temperature, for an attractive uni- 
tary Fermi gas (upper panel, Fig. 2a) or for a repulsive 
unitary Fermi gas (lower panel, Fig. 2b). The exact free 
energy would be bounded by the two successive expan- 
sion predictions and therefore could lie in between. With 
increasing the spin imbalance, the free energy of interact- 
ing Fermi gases approaches to the ideal, fully polarized 
limit. The reduction of interaction energy is clearly visi- 
ble. For the case with repulsive interactions at T = 5Tp, 
where we believe the quantitative applicability of quan- 
tum virial expansion, the minimum of the free energy 
occurs at the balanced limit of P = 0. This implies that 
there is no itinerant ferromagnetism at such a large tem- 
perature 5Tp. It therefore gives an upper bound for the 
ferromagnetic transition temperature of a homogeneous 
Fermi gas at unitarity. 



B. Free energy of a trapped unitary Fermi gas 

In analogy with the uniform free energy in Figs. 1 and 
2, we show in Figs. 3 and 4, respectively, the dependence 
of the free energy on temperature and on spin imbalance. 
Due to the existence of harmonic traps and hence a factor 
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Figure 2: (color online) Free energy of a homogeneous uni- 
tary Fermi gas as a function of spin imbalance at a given 
temperature (as indicated), for attractive interactions or for 
repulsive interactions. The exact free energy could lie between 
the second- and third-order virial expansion predictions. 



of rt~ 3 / 2 reduction for the n-th virial coefficient as shown 
in Eq. (|23[) . the applicability of quantum virial expansion 
is greatly extended. At the third order, we estimate that 
the virial expansion is now quantitatively reliable down 
to 0.5Tp for an attractive Fermi gas and down to Tp for 
a repulsive Fermi gas. 

We estimate also the upper bound for the ferromag- 
netic transition temperature of a trapped repulsively in- 
teracting Fermi gas at unitarity. As seen from Fig. 4b, 
the trapped free energy does not exhibit any signal for 
itinerant ferromagnetism at the Fermi degenerate tem- 
perature Tp. We thus conclude that (T c ) upper < Tp for 
a trapped unitary Fermi gas. 



Absence of itinerant ferromagnetism in the 
third order virial expansion 



The absence of itinerant ferromagnetism in our virial 
expansion may be analytically understood. To the third 
order, we calculate the two by two susceptibility matrix 
S = (dn a /dfi a >). A divergent or negative spin suscep- 
tibility xs = (dSn/dS/i) in the balanced limit of P = 
would signal the instability to the formation of ferromag- 
neitc order. The calculation of the susceptibility matrix 
and xs is straightforward. We find that for a homoge- 
neous Fermi gas at unitarity, 



S(P = 0) = 



1 



k B TX 3 



A B 
B A 



(30) 
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Figure 3: (color online) Temperature dependence of the free 
energy of a trapped unitary Fermi gas. We use the same no- 
tations as in Fig. 1. The experimental data for an attractive 
trapped Fermi gas at unitarity, reported by Salomon group 
at ENS [20l |. are shown by squares. 
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Figure 4: (color online) Free energy of a trapped unitary 
Fermi gas as a function of spin imbalance at a given tem- 
perature as indicated. The same notations as in Fig. 2 have 
been used. 



where 



A = —= 



Viz 



'it J (1 + ze-*) 
o 



rdt + 2z 2 Ab 2 , H + 5z 3 Ab 3M (31) 
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Figure 5: (color online) Spin susceptibility of a homogeneous 
Fermi gas at unitarity, normalized by the T = ideal Fermi 
gas susceptibility value xo = 3n/(2ep), with n the density 
and tF the Fermi energy. The spin susceptibility of a repul- 
sive or attractive Fermi gas at unitarity lies above or below 
the ideal Fermi Gas result (dot-dashed line). The inset shows 
the inverse of the spin susceptibility of a repulsive unitary 
Fermi gas. A linear extrapolation of the high-T virial expan- 
sion result at T/Tf — 4 ~ 6 to low temperatures indicates a 
ferromagnetic transition temperature of T c ~ Q.42V. 



6Tp) to low temperatures, we estimate very roughly that 
T c ~ 0AT F . 

The spin susceptibility of an attractive unitary Fermi 
gas shown in Fig. 5 is also of interest. At finite temper- 
atures, it is related to the measurement of the thermal 
spin flucutations: 

. fere. (34, 

N n 

A shot noise measurement of the spin fluctuations there- 
fore could be used as a sensitive thermometry for strongly 
interacting Fermi gases [35], provided that the spin sus- 
ceptibility is known. This seems to be now accessible, 
since the shot noise measurements of the density fluctua- 
tions in a weakly interacting Fermi gas have already been 
demonstrated very recently [361 |37| . We find that the 
spin susceptibility of an attractive Fermi gas at unitarity 
is strongly suppressed with respect to the ideal Fermi gas 
result, even well above the degenerate temperature Tp. 
At T = T F , the reduction is about 40%. 



and 



IV. CONCLUSIONS AND REMARKS 



B = 2z 2 Ab 2 H + 4z 3 Ab 3 ,H- 



(32) 



The spin susceptibility xs = 2(A - B)/(k B TX 3 ) is then 
given by, 



k B TX 3 



Viz 



'Jrj (l + ze-*) 
o 



z 3 Ab 



3,H 



(33) 



It is easy to see that, because of a positive third virial 
coefficient Ab^fj, the spin susceptibility of a repulsive 
unitary Fermi gas would always be finite and positive, 
for a finite fugacity. Therefore, there must be no fer- 
romagnetic transition if we restrict the quantum virial 
expansion up to the third order. The same observation 
is found for a trapped repulsive Fermi gas at unitarity. 



D. Third order expansion of the spin susceptibility 
for a homogeneous Fermi gas at unitarity 

Figure 5 reports the numerical result of Eq. (|33[) for a 
homogeneous repulsive and attractive Fermi gas at uni- 
tarity, where the fugacity z is solved consistently to the 
third order expansion in the number equation. 

The spin susceptibility of a repulsive unitary Fermi gas 
is larger than that of an ideal, non-interacting Fermi gas, 
as expected. Accordingly, its inverse Xo/xs is smaller 
(see the inset). This might be a high-temperature in- 
dication of the ferromagetic instability. The transition 
temperature is conveniently determined from the condi- 
tion Xs 1 (Tc) = 0. By linearly extraplating the high- 
temperature result of x^ 1 ( m the window of T — 4Tp ~ 



In conclusions, we have presented a quantum virial ex- 
pansion theory for the thermodynamics of strongly in- 
teracting, spin-population imbalanced Fermi gases. The 
(imbalanced) virial coefficients are calculable through the 
cluster partition function with different spin configura- 
tions. We have shown the virial expansion of thermody- 
namics up to the third order and have discussed in de- 
tail the numerical results for the high-temperature free 
energy of a homogeneous or trapped Fermi gas at uni- 
tarity. The predicted free energy may confront future 
high-precision thermodynamic measurements. 

An interesting application of our virial expansion the- 
ory is the determination of an upper bound for the fer- 
romagnetic transition temperature of a repulsive Fermi 
gas at unitarity. We have proven numerically and ana- 
lytically that, up to the third order quantum virial ex- 
pansion fails to predict itinerant ferromagnetism. As the 
third-order virial expansion is quantitatively applicable 
at T > Tp for a trapped repulsively interacting Fermi 
gas at unitarity, we therefore have estimated an upper 
bound (T c ) upper — Tp for the critical ferromagnetic tem- 
perature in harmonic traps. Further improvement of the 
upper bound requires the calculation of higher-order (im- 
balanced) virial coefficients, by using few-fermion solu- 
tions as the input (3^ - i40j . 

Much wider applications of the quantum virial expan- 
sion method for strongly interacting quantum gases seem 
to be feasible. In the near future, we are interested in 
inve stig ating the virial expansion of the Tan's contact 
[ill . l42j | the three-body recombination rate [43l |44| , and 
the quantum viscosity (4|| [46[ . 
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